Non-Zeeman splitting for a spin-resolved STM with a Kondo adatom in a 
spin-polarized two-dimensional electron gas 
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We theoretically investigate the spin-resolved local density of states (SR-LDOS) of a spin-polarized 
two-dimensional electron gas in the presence of a Kondo adatom and a STM probe. Using Green 
function formalism and the atomic approach in the limit of infinite Coulomb correlation, it is found 
an analytical SR-LDOS expression in the low temperature regime of the system. This formal result 
is given in terms of phase shifts originated by the adatom scattering and Fano interference. The 
SR-LDOS is investigated as a function of the probe position and different Fano factors. Our findings 
provide an alternative way to spin-split the Kondo resonance without the use of huge magnetic fields, 
typically necessary in adatom systems characterized by large Kondo temperatures. We observe a 
non-Zeeman spin-splitting of the Kondo resonance in the total LDOS, with one spin-component 
pinned around the host Fermi level. Interestingly, this result is in accordance to recent experimental 
data reported in Phys. Rev. B 82, 020406(R) (2010). 

PACS numbers: 72.25.-b, 73.23.-b, 74.55. +v, 75.20.Hr 



I. INTRODUCTION 

The scattering of electrons by a magnetic impurity in 
a metallic environment is responsible for the manifesta- 
tion of the Kondo effect X This phenomenon occurs as a 
result of an antifcrromagnctic coupling between the lo- 
calized spin at the impurity and the surrounding con- 
duction electrons of the host. In particular, at tem- 
peratures much lower than the Kondo temperature Tk, 
an electron cloud emerges to screen the magnetic mo- 
ment placed at the impurity site. Thus, a sharp reso- 
nance pinned at the Fermi energy appears in the impu- 
rity density of states and characterizes the formation of 
the Kondo peak. Such effect was first observed in resis- 
tivity measurements of magnetic alloys, later in trans- 
port properties of quantum dots (QDs) performed in a 
two-dimensional electron gas (2DEG). 2 ~— More recently, 
Kondo effect was also measured using scanning tunneling 
microscope (STM) in the presence of impurities deposited 
on metallic surfaces*^ - — 

In the context of unpolarizcd scanning tunneling mi- 
croscope (STM) probes, the conductance exhibits the 
Fano line shape due to the quantum interference between 
the transport channels given by the conduction bands 
and the adatom^i - — Additionally, for STM probes not 
very close to the metallic host and in the low tempera- 
ture regime, the STM device probes the local density of 
states (LDOS) of the sample. 

In the case of spin-polarized STM probes, interesting 
new features emerge as the spin-splitting of the Fano- 
Kondo profile of the conductance and the Fano-Kondo 
spin-filter i22r— In this scenario, several experimental and 
theoretical works discuss related phenomena employing 
ferromagnetic leads coupled to QDs and adatoms^ir— . 
In particular, in the emerging field of spintronics, the 



interplay between the Fano-Kondo effect and the ferro- 
magnetism of a metallic environment, plays a crucial role 
in the manufacturing of novel spintronics devices. 



Unpolarized STM probe 
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FIG. 1. (Color online) STM device composed by an unpolar- 
ized probe and a Kondo adatom hybridized to a SP-2DEG. 
The parameters V, tdR and w are hopping terms. The letters 
d a , *P,t and denote, respectively, the fermionic annihi- 

lation operators to the adatom, to the site of the host side- 
coupled to the adatom and to the site just below the probe. 

In this work, we report an analytical description of the 
spin- resolved local density of states (SR-LDOS), for an 
unpolarizcd STM probe, hybridized with a single Kondo 
adatom in a ferromagnetic (FM) island, considered here 
as a spin-polarized two-dimensional electron gas (SP- 
2DEG). Such analysis is done in the framework of the sin- 
gle impurity Anderson model (SIAM)££, using the atomic 
approach with infinite Coulomb correlation^— in order 
to determine the adatom Green function (GF). 

The main result of our simulations is the emergence of 
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an asymmetric spin-splitting in the Kondo peak, which 
agrees with some experimental worksi 34 ' 45 ' 47 ' 49 ' 50 Here 
we highlight the experiment done by S. L. Kawahara 
et ai 4 ^, where the measured LDOS shows that one spin 
channel has a Kondo peak pinned at the vicinity of the 
host Fermi level, while the opposite is shifted by the 
spin-polarization of the 2DEG forming the island surface. 
We note that the behavior we see here is not similar to 
the usual Zeeman splitting due to an external magnetic 
field^— — — i22i In the present model there is no external 
fields. 

Contrary to the QD systems, where the Kondo tem- 
perature is of the order of milliKclvin and the temper- 
ature could be tuned to observe the suppression of the 
Kondo resonance, in most STM systems,— the Kondo 
temperature is of the order of tens of Kelvin and the 
experimental STM setups do not allow the temperature 
variation in such range to observe the evolution of the 
Kondo effect. On the other hand, to verify the splitting 
of the Kondo peak in STM experiments, it is necessary 
to apply a magnetic field of hundreds of Tcsla, which is 
not feasible. Thus, the type of experiment proposed in 
our work could be useful to produce a splitting of the 
Kondo resonance without huge magnetic fields. 

This paper is organized as follows. In Sec. HH we 
show the theoretical model for the STM in terms of the 
Anderson Hamiltonian. We derive in Sec. IIII1 the SR- 
LDOS formula with and without a STM probe. For both 
cases in the low temperature limit, we show that the SR- 
LDOS expression can be labeled by phase shifts due to 
the adatom scattering and Fano interference. In Sec. IIV1 
we discuss the results of the host Fano parameter, which 
displays spin-polarized Friedel oscillations^ 1 - - — and we 
present a wide analysis of the SR-LDOS as a function of 
the bias-voltage, for different probe positions and Fano 
factors. The atomic approach is employed considering 
an infinite Coulomb correlation in order to calculate the 
adatom GF. We apply our formulation in Sec. |V] to de- 
scribe the Kondo peak splitting found in the experiment 
of Ref. [49|]. The conclusions appear in Sec. I VII and in 
the Appendices, we give details of the derivations of the 
atomic GF for the adatom as well as for the host Fano 
factor. 



II. THEORETICAL MODEL 

A. Hamiltonian 

In Fig. [I] we represent an unpolarizcd STM probe cou- 
pled to the FM island hybridized to the Kondo adatom 
deposited on its surface. Note that when the hopping 
term tdR >> w, the setup behaves as a Single Electron 
Transistor (SET^i 3 -, which we call peak limit, due to the 
emergence of the Kondo resonance in the LDOS energy 
profile as we shall see. The other limit we call dip limit, 
it comprises the cases tdR « w and tdR — W, which 
resemble the T-Shaped QD device 7 - characterized by a 



Fano-Kondo dip. 

The system we investigated is described according to 
the Hamiltonian 



H = H 



FM 



■H 



probe 



(1) 



The first term represents the SIAM 5 ^ 
Hfm = E £kljC \, 



ka k(7 



E 

ka 



Vdka(c\ d 
v ka 



d t c ka) + Ud\d^d\di, 



(2) 



that assumes the island as a SP-2DEG described by the 
operator ct (cj CT ) for the creation (annihilation) of an 
electron in a quantum state labeled by the wave vector 
k, spin a and energy 



£fc<7 — Do-kpQ. (k kp a 



(3) 



that depends on the spin-polarized band half-widths Do- 
forming the host Fermi sea and the wave numbers fcjr ff 
evaluated at the Fermi energy ep = 0. 

For the adatom, d\, (d a ) creates (annihilates) one elec- 
tron with spin a in state Ed- The third term hybridizes 
the adatom level and the host continuum of states. The 
coupling matrix element is modeled as 



Vdko 



V 



r 2 



FMt 



T 2 



(4) 



which obeys a Lorentzian behavior for a sake of simplic- 
ity to mimic a nonlocal coupling between the adatom and 
the island. NpMa is the number of conduction states for 
a given spin and T is the width of this interaction around 
£f- Coulomb correlation between two electrons with op- 
posite spins at the adatom site is also taken into account 
and is represented by the letter U . For simplicity we con- 
sider the U infinite version of the atomic approac h 57 ' 58 , 
that gives the adatom GF. In particular, taking the limit 
r 3> Ska in Eq. @, we obtain a constant Vdka, i.e., 



Vdka = 



This corresponds to the case of a site 



of the FM island side-coupled to the adatom, which we 
designate local coupling. We mention that Eqs. (J3j) and 
(U) were previously applied in the context of unpolar- 
ized bulk electrons in a system described by the Kondo 
model.— In this work, we employ the Lorentzian shape 
to emulate the nonlocality of the adatom-island coupling. 

The FM island is considered a spin-polarized electron 
bath, with polarization given by 



P 



pFMf + PFMi 



where 



PFMo 



2D n 



(5) 



(6) 



is the host density of states in the flat band approxima- 
tion. This quantity is expressed in terms of 



Do 



(1 + aPy 



(7) 
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for a given spin er and the unpolarized density 



written in terms of the width Do- The unpolarized part 
of the system is the conduction band given by the Hamil- 
tonian 

7~Lprobe ^ y ^pbfiybjia) (9) 

p*<7 

which corresponds to the second term in Eq. ([T]) for free 
electrons in the STM probe. Such conduction electrons 
are ruled by fermionic operators &t CT and bp a . To perform 
the coupling between Eqs. ([5J and (|9|), we have to define 

Kun = b k (V^R* + tdRda + H.C.) , (10) 
P 

as the spin tunneling Hamiltonian that hybridizes the 
STM probe conduction states with those in the island and 
the adatom. The former hopping parameter, we consider 
proportional to the energy independent term w, with the 
fermionic operator 

** CT = 7 ^5>^c fcV (11) 

which describes a conduction state at the site R later- 
ally displaced from the adatom. This operator admits an 
expansion in plane waves, due to the assumption of an 
infinite 2DEG forming the island surface. The second hy- 
bridization parameter in Eq. (|10[) is proportional to the 
adatom operator d a and the spatial dependent hopping 

t dR = t do exp(-k F R), (12) 

which provides a decreasing STM probc-adatom cou- 
pling. This characteristic ensures a vanishing behavior 
for huge lateral displacements, which was already em- 
ployed in the literatur e 14 ' 31 ' 32 . 

B. The atomic approach 

In order to implement the atomic approach^—, we 
consider the adatom-island coupling as local. Thus, we 
begin 

H F M = Efc CT £kac\ a C ia + Y, a E dX d ,ca 

+ V (<o^ + HXd,o a ) , (13) 

which is derived from from Eq. @ taking into account 
Vdka = ,*Y ■ Here X a h is the Hubbard operato r 64 ' 65 
that projects out the doubly occupied state from the 
adatom to ensure the limit of infinite Coulomb corre- 
lation, where the label (a, b) defines the parameters asso- 
ciated with the corresponding atomic transition and ty a 



is [Eq. (fTTj) ] evaluated at the origin. This Hamilto- 
nian is useful for the calculation of the adatom GF and 
consequently the system SR-LDOS. 

This GF is based on an extension of the Hubbard cu- 
mulant expansion also applicable to the Anderson lattice 
with impurity-host couplings treated as perturbations^ 
The use of this expansion allows to express the exact GF 
in terms of an unknown effective cumulant. In a pre- 
vious work, it was studied the Anderson lattice with an 
approximate effective cumulant obtained from the atomic 
limit of the model in a procedure called as the zero band 
width (ZBW) approximation. The advantage is that such 
method includes all the higher order cumulants absent in 
the previous diagrammatic calculations^ 5 - 

The method presents some similarities to the exact di- 
agonalization (ED) for the SIAM. The ED is a brute- 
force method to solve the Hamiltonian treated as a dis- 
crete bath by considering the impurity coupled to a fi- 
nite number of conduction sites (N c ) of the host band. 
In principle, it is an exact method as the name implies, 
but its limitation relics in the number of conduction sites 
considered, the Hilbert space grows extremely fast when 
N c is enlarged. The cumulant approach is not equivalent 
to the ED, it allows to derive the adatom GF from the 
Anderson lattice model in the atomic limit ^ This GF is 
related to the conduction band GF. which depends on the 
atomic level that represents the band. Such level is the 
starting point for the ZBW approximation. The ZBW 
is a mapping of the system Hamiltonian onto an effec- 
tive tight-binding chain with only two sites, one for the 
adatom discrete energy and other for the band atomic 
level. Then employing the Lehmann representation, the 
atomic GF for the impurity is calculated as well as the 
approximated cumulants. They determine the full GF in 
such a way that the Friedel sum rule should be satisfied 
by the band atomic level. In particular, this GF does 
not present a spurious oscillatory behavior in the Kondo 
peak as the observed by the standard ED method. In the 
ED context, this artifact can be solved by the self-energy 
trick^ We also would like to mention, that the atomic 
approach could be easily generalized to multi-orbital An- 
derson model, because once the atomic solution is known, 
the method follows the same steps realized in the spin 
5 = 1/2 SIAM. 

To obtain the exact GF of an Anderson impurity 
(adatom), we can employ the chain approximation^ 5 - but 
considering all the possible cumulants in the expansion 
for the Anderson lattice. Similarly to the Feynman dia- 
grams, it is possible to rearrange all those that contribute 
to the exact adatom GF by defining an effective cumu- 
lant, determined by all the diagrams that cannot be sep- 
arated by cutting a single edge ( "proper" or "irreducible" 
diagrams). 

As we are interested in the exact GF for the adatom, 
we use the standard definition 

Qf (r) = - l -e (r) Tr {g FM [d a (r) , 4 (0)] + } (14) 

in time coordinate, where h is the Planck constant h di- 
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vided by 27r, 9 (r) the step function at the instant t, Tr 
the trace over the eigenstates of the Hamiltonian in Eq. 
©; Qfm the density matrix of the FM island and [, ], 
is the anticommutator between the adatom operators at 
different times. 

The time Fourier transformation of Eq. ([M]) thus pro- 
vides the adatom GF in energy coordinate, which is then 
obtained by replacing the bare cumulant for the effective 
one calculated by following the atomic approach with the 
Hamiltonian in Eq. (fT3j) . As a result, we have 



l-M e J f {u)\V\*Y.%G C a{k,w) 



(15) 



as the adatom GF in terms of the effective cumulant 
M e Jf{ui) and the free-electron GF 



UJ ' £ ka + iT l 



(16) 



where rj — > + . The atomic version of Eq. (|15|) is given 
by (see Appendix A) 



CM 



i-M?(u)\v\*gz BW (u) ' 



which results in 



Git A*) 



i+GtiM\v?G% Bw {uy 



(17) 



(18) 



for the effective cumulant determined from the adatom 
GF calculated in Appendix [S] both dependent on 



■>ZBW 



i 



(e a - A 4 ) + V 



(19) 



for an electron state, in the ZBW approximation with 
fj, to denote the FM island chemical potential. As we 
can see, Eq. ([T^]) replaces all energy contributions of 
the original Fermi sea by two spin dependent atomic lev- 
els, i.e., we perform the substitution Y% e kij c\ cr —> 

— Ka KG Ka 

J2 a £ 0crcl a c 0a in Eq. ([131) witn £ ka = £oa representing 
the band atomic level for a given spin a. As this pro- 
cedure overestimates the coupling of the spin-polarized 
conduction bands of the island with the adatom due to 
the concentration of the bands at atomic levels, we have 
to moderate this effect^ performing the substitution of 
V 2 by A 2 in Eqs. (JT7J) and (fig]), where A = ttV 2 p is 
the Anderson parameter. 

To determine the adatom GF, we use the atomic cu- 
mulant Ad^ioS) as effective in Eq. ([15)) and verify that 



Gi d (uj) 



1 



(20) 



provides an analytical expression in the flat band approx- 
imation. We mention that .M" (a;) is a simplification, 
but it contains all the diagrams that should be presented 



in such way that the correspondent GF displays realistic 
features. 

As the final step of the atomic approach implementa- 
tion, we have to find adequate values of the atomic levels 
e 0cr that well describe the ZBW GFs in Eq. (HI]) and 
consequently the adatom GF. To that end, we use the 
condition that, in metallic systems, the most important 
region for conduction electrons is placed at the Fermi en- 
ergy ef and that the coupling to an Anderson impurity 
leads to the Friedel's sum rule^i 



Pd,a( £ F) 



1 



{Gi d { £ F)} 



sin 2 (8 a {e F )) 
7rA„ 



(21) 



for the adatom spectral density. Here S a (sF) = "Kn^a is 
the conduction phase shift at the Fermi level, 3 repre- 
sents the imaginary part, A CT = A(l + aP) is the spin 
dependent Anderson parameter and Ud,a is the adatom 
occupation with spin a. Thus we find the atomic levels 
Equ calculating self-consistently Eq. ((2"Tj) together with 



nd,cr =< Xd.uu >— 

7T 



duS{Gi d (u))}n F (uj). 



(22) 

In Eq. (f2"2"|) . rip (to) is the Fermi-Dirac distribution. 

It is important to emphasize here that the choice of 
the atomic approach to calculate the adatom GF is only 
due to its computational simplicity and ability to obtain 
the Kondo peak, but we must take into account that the 
method presents some limitations that were extensively 
discussed in the original papers^r— . The SR-LDOS for- 
mulas obtained in Sec. IIIII are general, we could employ 
others more powerful and precise techniques to calcu- 
late the GF of the Anderson impurity, like the Numer- 
ical Renormalization Group (NRG ) 10 ' 11 ' 68 without any 
modifications in the formalism. 



III. SR-LDOS 



A. SR-LDOS in the scheme of phase shifts for the 
FM island 



In this section we derive at the temperature range 
T <C Tk, the SR-LDOS for the FM island with an 
adatom following the procedure found in Ref. [62| , which 
was applied in the Kondo model with unpolarized bulk 
electrons. Such method allows to express the SR-LDOS 
in terms of the phase shifts due to the adatom scattering 
and the Fano effect. This latter is originated in the in- 
terference between the electron paths formed by the host 
conduction band and the Anderson impurity. We empha- 
size that the STM probe is not present in this section. 
Initially we derive a formalism for finite T and later on 
we take the limit T — > 0. It is well known that 



Pldos ( w > r ) 



1 



(23) 
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provides the SR-LDOS formula. The GF g a (uj,R) is 
obtained from the Fourier transform of 

G a (r, R) = ~6 (r) Tr j g FM [«£ Ra (r) , tf Ra (0) j ^ 

JH (24) 

in time coordinate, where Qfm and [, ] , are the density 
matrix of the FM island Hamiltonian and the anticom- 
mutator between the operators given by Eq. ([lip at dif- 
ferent times, respectively. 

In order to be explicit, we have to apply the equation 
of motion procedure (EOM) on Eq. (|24| to demonstrate 
that Q a (u),R) is coupled to other GFs as follows 

G a (u>,R) = g„ (w, 0) + g a (u, R) % (w) g a (w,-R). (25) 

The first term 



g a (w, R) 



1 



E — 



(26) 



describes the bare GF for an uncorrelated electron state 
at the site i?, laterally displaced from the adatom and 



9a (w,i?) 



1 



E 



(27) 



is the correspondent one dressed by the nonlocal hy- 
bridization [Eq. (QJ]. As a scattering center, the adatom 
defines a scattering amplitude 



%(<*>) = — g?(w) 

npo 



(28) 



proportional to the GF Q'jf (uj). 

The emergence of Fano interference and Friedel oscilla- 
tions in the system SR-LDOS is a result of the interplay 
between Eqs. (f2l)| and (f2"T)) . These effects can be eluci- 
dated by regrouping the terms in Eq. (|25|) in such way 
to achieve the form 



deduced from Eqs. © and ©. The SR-LDOS also de- 
pends on the function 



A a (R) = —%{g a (u,R)} 
npo 



PFMa T n u , T 2 
-JO {KFaR) ■ 



r 2 



-UJ Z 

(32) 

that encloses spin-polarized Friedel oscillations for 
charge. Details of the method employed to derive Eq. 
pop appear in Appendix |B] We mention that to compare 
the SR-LDOS simulations of Sec. |V]to the experimental 
data measured in a Fc island with a Co adatom^, it is 
necessary to set the Fano parameter in Eq. (|30|) . as a 
spinless ratio and equal to zero. 

At low temperature limit, the SR-LDOS is well de- 
scribed by Eq. (|29|) at the ground state. In such region, 
it is possible to classify the SR-LDOS in terms of the spin- 
dependent phase shift 8 a (oS) for the conduction electrons 
due to the adatom scattering center. According to the 
results on the SIAM§£, this characterization is obtained 
using the relation 



exp[2i8 a (u)\ = 1 - 27:p ow ^-J2k ( pq? ) 7 ^( UJ ) i 



X 8 (U) — Ska) 



(33) 



that correlates the phase shift 8 a (w) to the real and imag- 
inary parts of the scattering amplitude T a (w), given by 
Eq. (|28p. thus resulting in the formula 



tan (So- (w) 



3{T ff M} 



(34) 



For the Fano interference we define an analogous 
relation,— ~— introducing the phase shift 8 qFM as 



, c ^{9a (U,R)} qFMa 



in order to show that 



%{g* (",R)} : 



(35) 



Pldos ( w > R) = PFMa + ttpI { (A 2 a (R) - q FMa ) 



PLDOS ( W > R) = PFMa [l - Jo {kFaR) Fa (w)] , (36) 



x 3 {% («)} + 2A a (R) qpMaVi {% (w)}} , (29) 
where 5i means real part, 

qFMa = 3? {g a (ui, R)} = PFMrT J (k F aR) ^r— — 

■Kpo po T z + uj z 

(30) 

represents the Fano parameter due to the adatom-island 
hybridization in the wide-band limit and characterized 
by spin-polarized Friedel oscillations in the zeroth-order 
Bessel function J (kp cr R). The spin-dependent Fermi 
wave numbers kp^ and kp-f are related to each other via 



\/ TTp kF ^ 



(31) 



becomes the expression for the system SR-LDOS and 
characterized by 

t- / s , cos 2 (S a (u)-8 qFM ) 

^ M = 1 ' (37) 

exclusively expressed in terms of the adatom scattering 
and Fano phase shifts. The last formula is the main result 
of this section, it represents the SR-LDOS of a metallic 
surface considered as a spin-polarized 2DEG coupled via 
a nonlocal hybridization to an adatom, in the framework 
of the SIAM given by Eq. ©. In Sec. ITlTBl we shall 
see that the STM spin-resolved conductance formula, as- 
sumes identical structures to Eqs. (|3^|) and (|3T|) . but with 
a redefined Fano parameter due to the probe presence. 
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B. Differential conductance and the SR-LDOS 
probed by the STM 

In the low temperature regime and for a STM probe 
not very close to the metallic sample, the spin-resolved 
and differential conductance of the STM device is pro- 
portional to an effective SR-LDOS. 

To derive this effective SR-LDOS, we have to imple- 
ment the linear response theory treating the tunneling 
Hamiltonian in Eq. (|10p as a perturbation, just to en- 
sure the weak tunneling regime as verified in experimen- 
tal conditions^. Taking these assumptions into account 
we show that the spin-current follows the expression 

2e f 

I a = j-^w / duj [n F (w - e<f>) - n F M] p LDO s M i 

. (38) 

with T w = 2tt \ w\ pQ as the parameter that hybridizes 
the system conduction bands, e the electron charge, <f> 



the bias- voltage and 



Plnos (u,R) = ~$t{g„ (w,iZ)} 



(39) 



as the effective SR-LDOS probed by the STM device. 
Such density is calculated using the GF Q a (u, R) ob- 
tained from the time Fourier transform of 



(r),*L (0) 



Q a (t,R) =-jj0(r) Tr\ QFM 



written in terms of the operator 

^Ra = ^Ra + {vAp ) l/2 q R d a 



(40) 



(41) 



that describes the couplings between the probe and the 
island with the adatom, which depends on the new Fano 
parameter defined by 



<1R 



{irApo) 1/2 (t dR /w) = q cxp(~k F R), (42) 



due to the interference between these additional conduc- 
tion channels. To obtain the differential conductance 
G a = -§^I<r for a given spin we consider Eq. (|38|) and 
show that 



2e 2 

G a = —, — 
h 



duj 



d I 



<P)j Pldos (u,R) 
(43) 

is the spin-resolved conductance for the STM device. 

Note that the spin effects on the system conductance 
lie on the free density of states of the island in Eq. © 
and on the scattering amplitude T a (ui) in Eq. (|2"5)l due 
to the adatom. Thus, we need to express Eq. (|39p in 
terms of the adatom GF g dd (w) by employing the EOM 
procedure. This method leads to 

Q a (u, R) = Ga (w, R) + 7rAp 4^ d (u) + (7rApo) 1/2 qn 
x g* d (cj, R) + (nApof 2 q R Gf (w, R) , (44) 



and displays that such GF is linked to the Eq. (l25t for 
the site R of the island and simultaneously to 

g* d («, r) = (ttA Pq ) 1/2 [ qFMa - iA a (R)} gf («) , 

(45) 

and g^ (uj,R), where the former expression is obtained 
from the time Fourier transformation of 

g^ d (r) = - l -e (r) Tr {g FM (r) , 4 (0)] + } , (46) 

which is also equal to g^ d (uj,R). Thus replacing Eqs. 
|, gll) and 



f^Po\ 



1/2 



{A a {R)^{gi d h} 



\ 7T / 

- q FM ^{g d a d {L0)}}, (47) 
into Eq. ([3U]). we find that 

Pldos ( w > r ) = PFMo + Trpo { (Al (R) - q Ra ) 

x Ss {% (w)} + 2A a (R) q Ra $t {% (w)}} 

(48) 

obeys Eq. (f2U)) for the SR-LDOS in the case of an absent 
STM probe, but with an effective Fano parameter 



qRcr = qFMa + 



(49) 



that takes into account an intrinsic Fano interference 
(qFMa) due to the adatom-island coupling represented 
by the first term and an extrinsic one (q R ) as a result 
of the STM probe hybridized with both adatom and the 
metallic surface, enclosed by the second part. We pointed 
out that for T <C Tjc, the spin-resolved conductance in 
Eq. (|43)) becomes directly proportional to the effective 
SR-LDOS evaluated at the energy e</>. To show that, we 
use the Dirac delta distribution expressed by the minus 
derivative of the Fermi function in Eq. (|43[) , which elim- 
inates the integration over energy and gives 



Ga 



2e : 



K^wPlDOS ( e <? 



(50) 



This result means that the effective SR-LDOS is a fairly 
representative function for the spin-resolved conductance 
of the system, which favors us to apply the zero temper- 
ature formalism of phase shifts discussed in Sec. IIII Al bv 
introducing 



, j- qria 
tano„„ = — - — ; — -, 
qRa A a (R) ' 



(51) 



as the total Fano phase shift. Combining it with Eq. 
(p4| for the scattering amplitude T a (w), it is possible to 
derive 



Pldos ( w > r ) = Pfmo 1 - Jo ( k FoR) Fa (w, R) 
and 



Ta (w, R) = l- 



cos 2 (8„ (w) -S qR J 



cos 2 6, 



(52) 



(53) 



qRa 
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as expressions that represent the effective SR-LDOS 
probed by the STM. The latter contains two sources for 
Fano effect: the first concerns on the Fano interference 
between the traveling electrons through the host conduc- 
tion band that can "visit" the adatom site and go back 
to it, and those that do not perform such "visit". Addi- 
tionally, the second process is composed by the couplings 
of the probe with the island and the adatom. It also has 
the scattering of the traveling electrons due to the side- 
coupled adatom, which in certain conditions leads to the 
Kondo effect. We also remark that the phase shift 5 a (10) 
given by Eq. Q34p . is obtained in this work employing the 
atomic approach. 

From Eq. (f52"j). we are able to determine the following 
occupation number 

/+oo 
dupLDOsfa, R = 0)n F (w). (54) 
-oo 

As we shall see, this formula will guide us to better un- 
derstand the results of Sec. IIVI 



IV. RESULTS 

Here we present the results obtained via the formu- 
lation developed in the previous section. We employ a 
set of parameters that defines the Kondo regime: Ed = 
-10.0A and V = 12.0A, where A = 0.01£> (A) = 1-0). 
We will basically compare two regimes of qo value. The 
large qo limit we call peak limit due to the formation of 
a Kondo peak in the LDOS. The other regime, corre- 
sponding to small and intermediate values of qo, we call 
dip limit since the LDOS presents a dip around the Fermi 
level. 

It is necessary to mention that the SR-LDOS formula 
derived in the previous section is valid for zero temper- 
ature. On the other hand the numerical procedure of 
the atomic approach is well established for T -^i Tk but 
still finite T. So we consider a very low temperature 
T = 0.001A to allow the combination of both proce- 
dures. To investigate the spatial behaviors of the Fano 
factor and the LDOS, we define a dimensionless param- 
eter fcp^i? = kpR, to represent the STM probe-adatom 
lateral distance in Fig. [T] 



A. Intrinsic Fano parameter 

Here we look with some more detail the spatial de- 
pendence of the intrinsic Fano factor qF Kia • Fig. HJa) 
shows qpua against kpR for P = 0.3, ui = 0.02D and 
r = O.lui. Interestingly, this Fano factor reveals spin- 
polarized Friedel oscillations. These oscillations exhibit 
enhanced amplitudes for the spin-up channel and a phase 
shifted pattern in relation to the spin-down component, 
due to the spin-dependent Fermi wave number [Eq. (|31[) ] . 
This phase shift yields irregular oscillations in the total 
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FIG. 2. (a) Spin-polarized Fano factors qFMa for T = 0.1 
and (b) the sum qfm = qFMf + qFMi for T = 10.0, F = 20.0 
and r = 40.0 in units of to as a function of the dimensionless 
parameter /cf-R. 



Fano parameter qFM = QFMt + qFMi [Figs. Ela) and 
lib)]. 

In Fig. EJb) we can also observe by changing from 
r = 10cj to r = 40w, that the spin-polarized Fano pa- 
rameters approach to zero, exhibiting flattened oscilla- 
tions as a function of kpR- For large enough T, we have 
qpMo- ~ 0. Since the present atomic approach is valid 
only for constant Vdka which is obtained for large T (local 
coupling) , we will restrict our analysis to negligible qFMa 
values. In particular, for a perfect local coupling config- 
uration, settled by the condition qFMa — in Eq. (|30|) . 
we can conclude according to Eq. (|4U)) . that the Fano 
parameter q a induced by the STM probe is the only one 
that rules the Fano interference in the system. In Sees. 
IIV Bl and IIV CI we discuss the possible interference lim- 
its for q a in the local coupling regime, where the atomic 
approach is applicable. 
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B. SR-LDOS in the peak limit (q a = 100) 



Figs. Gfl»-(c) show the SR-LDOS and total LDOS for 
increasing polarization degree of the host. In the unpo- 
larized case P = (not shown) we have the well known 
twofold degenerate Kondo peak as in the SET^^ How- 
ever, for increasing P this degeneracy is broken and the 
Kondo peak splits into two peaks. While the correspond- 
ing spin- up peak presents a slight shift for increasing P, 
the spin-down one remains pinned around the host Fermi 
level. Additionally, the widths of these peaks show op- 
posite behaviors. While the up-peak broadens the down- 
peak shrinks as P enlarges. In the insets we present with 
more detail the SR-LDOS around the Fermi level, where 
we can see more clearly the spin-splitting of the Kondo 
resonance. 

Thus, it is valid to mention that our spin-splitting re- 
sembles the work of Y. Qi et al.^ for a spin current in- 
jected in a nonmagnetic conductor with a Kondo adatom. 
In this work, it is possible to note a tendency of a pinning 
for one spin component of the Kondo peak. Kondo peak 
splitting was also observed in a QD system hybridized 
to ferromagnetic reservoirs^ and in the presence of spin- 
flrp25. 

In the experimental point of view the spin-splitting of 
the Kondo resonance has already been observed in sys- 
tems of QDs coupled to ferromagnetic reservoirs^ and 
in a carbon nanotube QD interacting with a magnetic 
particle.— Observe that for P = 0.2 and P = 0.5 it is not 
possible to resolve the spin-splitting of the Kondo peak 
in the total LDOS. Experimentally, non-resolved Kondo 
peak splitting can also occur.— Although the STM ex- 
periments are in general restricted to energies around the 
host Fermi level, here we show the SR-LDOS for a wider 
energy window in order to see the polarization effects on 
the adatom level Ed, as displayed in Fig. [3] 

The analysis of the non-Zeeman splitting in the full 
LDOS for different STM probe positions is presented in 
Fig. H for P = 0.5. We see that the LDOS profile be- 
comes flatter for large enough values of the dimensionless 
parameter kpR, thus revealing a crossover from the peak 
limit at R = to the background value represented by 
the host free density of states [Eq. ((§)]. 

We end this section presenting Fig. [5J where the de- 
pendence of the spin-splitting of the Kondo peak AE as 
a function of the host polarization P can be observed. 
Note that AE displays a nonlinear behavior as P in- 
creases. This nonlinearity was also found by Y. Qi et 
alM 






FIG. 3. Spin-resolved local density of states (SR-LDOS) at 
R = 0, in arbitrary units (a.u.), as a function of the energy 
uj for increasing P value: (a) P = 0.2, (b) P — 0.5 and (c) 
P = 0.8. In the insets we show the SR-LDOS around the host 
Fermi level, where we can see the Kondo peak splitting with 
the pinning of the spin down component. We also observe 
that as P increases the spin down resonance shrinks. 



C. SR-LDOS in the dip limit (q < 1) 

In contrast to the qo = 100 case previously analyzed 
in the peak limit, here we discuss the qo < 1 regime. 
This corresponds to a direct STM probe and host sur- 
face tunneling being dominant. In this situation, the sys- 
tem reduces to an adatom side-coupled to the FM island 



as represented in Fig. [TJ which is equivalent to the T- 
Shaped QD device^. For small Fano factors (qo = 0.01), 
the LDOS exhibits antiresonances (dips) instead of reso- 
nances, as we can see in Fig. |6ja). This happens in such 
a way that the antiresonance in the up SR-LDOS chan- 
nel, disappears for high enough polarization as seen in 
the inset of Fig. [H^a). On the other hand, the dip of the 
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FIG. 4. Local density of states (LDOS), in arbitrary units 
(a.u.), as a function of the energy u for different UfR values. 
In the inset, we present the LDOS in detail. 
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FIG. 6. Spin-resolved local density of states (SR-LDOS) at 
R = 0, in arbitrary units (a.u.), as a function of the energy u. 
In the full curve we show the local density of states (LDOS). 
In the insets, we present the SR-LDOS and the LDOS in the 
vicinity of the Fermi energy oj ~ ef =0. For P = 0.2, there is 
a Fano-Kondo structure in the LDOS profile, but for P — 0.8 
an enhanced Kondo peak appears in the spin-down channel. 



FIG. 5. Kondo peak splitting AE in units of -Do, as a function 
of the host polarization P. 



down component shrinks as P increases. The width of 
this antircsonancc lowers two orders of magnitude when 
we change the polarization from P = 0.2 to P = 0.8. So 
we can conclude that the polarization induces a contin- 
uous second order insulator-metal transition in the sys- 
tem. As the polarization grows from P = to P = 0.8, 
the up SR-LDOS component becomes continuously flat, 
generating a finite SR-LDOS in the vicinity of the Fermi 
energy u ~ Ef = 0. At the same time, the down SR- 
LDOS component shrinks, closing the total LDOS gap, 
as indicated in the inset of Fig. H^a) for P = 0.8. 

The last cases in the dip limit we investigate, corre- 
spond to q a = 1, P = 0.2 and P = 0.8 as presented in Fig. 



[6^b). For the low polarization P = 0.2, the Fano inter- 
ference is robust presenting the coexistence of a structure 
composed by dips and peaks (Fano line shape). However, 
for the large polarization P = 0.8, the Kondo peak in the 
adatom GF prevails in the spin-down channel and sup- 
press the destructive interference. This feature appears 
in the insets of Fig. EJb). 

The simulations presented resemble experimental re- 
sults found in ferromagnetic contacts as in the M. R. 
Calvo et al. work^i (see their Fig. lc) and adatoms sys- 
tems discussed by N. Nccl et al£^ (see their Fig. lb). 
In the former experiment, depending on the electrodes 
composition coupled to the contact, the energy profile of 
the conductance displays patterns without resolved spin- 
splitting. This behavior is related to the cases of small 
and intermediate values of the Fano parameter, which 
were obtained with Fe and Co contacts respectively. For 
the second experiment, spin-polarized STM probes made 
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by Fe or W were employed on a Co adatom deposited on 
a Cu(lll) surface. In particular, in this case no resolved 
spin-splitting was verified, which corresponds to the limit 
of small Fano factor. These setups can be reproduced in 
our simulations considering intermediate polarizations. 



D. Phase shifts and occupation numbers 



In Sec. irriBl we show that the SR-LDOS can 
be expressed in terms of the phase shifts 5 a (to) and 
8 qiitT . Here we explore the quantities cos 2 5 qRa and 
cos 2 (5 a (cj) — 5 qRtl ) in order to gain further insight about 
the spin-splitting and Kondo peak pinning found in 
the SR-LDOS. In Fig. ([7]) we present these two func- 
tions against kpR for both spin components, where 
the numerical parameters arc indicated in the plots. 
The main distinction between the up and down curves 
can be observed in the range < fcp-i? < 5.0 where 
cos 2 (5 a (ef) — S qR<T ) w 1 for spin-down while it stays be- 



low 1 for spin-up. In particular, for kpR 
cos 2 (S qRtT ) — > which means that 5, 



'I n a 



— > we find 
7r/2. Conse- 
quently, cos 2 {5 a (ef) — 5 qRa ) — > sin z (5 a (ef))- On the 
other hand, according to Friedel sum rule [Eq. (|2"Tj) ] we 
have 5 a (ef) = Trader- This implies that rid,\, — > 0.5 while 
rid,-\ stays below 0.5. The fact that rid,\ remains close 
to 0.5 results in the pinning of the spin-down component 
of the Kondo resonance. For increasing kpR the curves 
obtained from Eq. ([BT]) display a series of peaks and dips 
due to the interplay between the exponential decay in 
Eq. (fT2|) and the Friedel oscillations. 

To complement the analysis on the emergence of the 
Kondo peak pinning, we present Fig. [5] In such fig- 



ure, we plot the occupation number 



LDOS 



ofEq. jH 



as a function of the Fano factor qo for differing spin- 
polarizations P. This gives also the same tendency of 

and rid a should pre- 



l LDOS 



nd,t an d nd,i, since both 

serve some proportionality. In the large qo limit, ^tuos 
approaches to 0.45 while n^ LDOS moves to lower values as 
P increases. This corroborates the pinning of the spin- 
down Kondo resonance. Note that the Friedel sum rule 
given by Eq. ([2"Tjl . ensures the pinning of the Kondo 
peak in the down-channel and the displacement of the 
up peak, via the inequality n^-j- < n>d,i ~ 0.5. In the 
inset of Fig. [5J the plots of rid,^ and na,i against P 
confirm this inequality. Such occupations indicate that 
there is a net magnetic moment at the adatom partially 
screened by the conduction electrons, where the lack of 
spin-down conduction electrons is not able to blind the 
spin-up component of the impurity, thus avoiding a for- 
mation of a defined Kondo peak in that channel. On the 
other hand, the excess of spin-up conduction electrons 
yields to a Kondo peak in the down-component of the 
adatom DOS. We can conclude that, there is no abrupt 
break down of the Kondo effect, but a crossover from the 
ordinary Kondo effect to a situation where the Kondo 
effect is gradually suppressed. 
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FIG. 7. Analysis of the system phase shift structure evaluated 
at the host Fermi level. Plots of cos 2 (<5 CT (ef) — S qRa ) and 
cos 2 8 qRtT , as a function of the dimensionless parameter kpR 
for both spin components. In the inset, we show the resonance 
and antiresonance features at kpR — 15.0. 



V. ANALYSIS OF SOME EXPERIMENTAL 
RESULTS 

In this section we successfully apply the present devel- 
oped formulation to reproduce recent experimental find- 
ings on a single Kondo adatom coupled to a magnetic 
cluster. To our best knowledge, the first experimental 
work that explores Kondo adatom on a ferromagnetic 
host was recently done by S. L. Kawahara et al^-, which 
used an unpolarized STM probe on top of a Fe island 
with a Co adatom. They observed that, depending on 
the adsorption site of the Co atom, a spin-splitting of 
a Fano-Kondo dip is induced by the spin-polarization of 
the island, and can be explained by a double Fano an- 
tiresonance in a single particle picture. 

As pointed out in Ref. [49[ there are two competing 
mechanisms that can result or not in the Kondo effect. 
The first one is the antifcrromagnctic exchange interac- 
tion between the Co adatom and the itinerant sp island 
electrons. The second one is the exchange interaction 
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FIG. 8. Occupation number rf LDOS of Eq. (|54|l as a function 
of the Fano parameter q and with different values for the 
spin-polarization P of the host. In the inset, we present n<i, CT 
against P. 



all the relevant parameter range of the SIAM in order to 
better reproduce the line shape of Kawahara et al. work. 
The optimized values are Ed = —3. OA and go = 0.01, 
which characterize the intermediate valence regime of the 
system. 
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due to the direct d — d ferromagnetic interaction of the 
Co adatom and the Fe atoms of the magnetic substrate. 
Our formulation only can be applied to this system if 
the antifcrromagnctic interaction dominates over the di- 
rect ferromagnetic correlation. This competing processes 
appear in Ref. [49| as an assumption, but further exper- 
imental and theoretical investigations are necessary to 
better understand the dominant mechanism. 

There are some related experimental results, supported 
by first principles calculations of M. R. Calvo et al£L, in 
related systems. They found the existence of a Kondo 
peak in ferromagnetic atomic contacts hybridized with 
electrodes (both built by Fe, Ni and Co), differently 
from those found in the bulk limit. In nanoscale, the 
electrons at the surfaces of these junctions experience 
interactions where the antiferromagnetic coupling over- 
comes the ferromagnetic correlations. In such setups, 
the nanocontacts play the role of the Co adatom used 
in the Kawahara et al work. Thus, our model in its 
present form does not support the opposite case charac- 
terized by strong ferromagnetic correlations, which are 
usually modeled by a Heisenberg type interaction. For 
an enhanced antiferromagnetic coupling between the sp 
electrons and the adatom, the picture of a spin-polarized 
electrons gas as discussed by M. R. Calvo et a/4£ can 
be employed to describe a ferromagnetic metallic sample 
with a Kondo impurity. Additionally, we would like to 
remark that Kondo adatoms and some QD systems in- 
deed have a spin S > 1/2, which can be detectable by 
a magnetic anisotropy signature 1 23 i 53 i 55 , I n these cases, 
a multi-orbital Anderson Hamiltonian could offer a more 
detailed modelin g 22 ! 51 and improve the accuracy of the 
present work. 

However, we changed the localized adatom level Ed in 
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FIG. 9. Spin-resolved local density of states (SR-LDOS) at 
R = 0, in arbitrary units (a.u.), as a function of the energy 
oj. In the full curve we represent the total density of states 
(LDOS). (a) In the inset (I), we present the splitting of the 
Fano-Kondo dip in the Kondo peak region for Ed — — 10A, 
whereas in the inset (II) we perform the same analysis for 
the intermediate valence case, using Ed = — 3A. In the inset 
(III), we present the suppression of the double Fano-Kondo 
dip structure as a function of kpR- 

It is clear from Fig. [HI that the resolved LDOS double 
dip structure is originated from the splitting of the up 
and down spin components. To show the evolution of the 
double structure at oj ~ £f — from the Kondo peak to 
the intermediate valence regime, we present in the inset 
(I) of Fig. © the splitting of the Fano-Kondo dip for 
Ed = —10 A whereas in the inset (II) we present the 
correspondent case in the same energy range, but for the 
intermediate valence situation, with Ed = — 3A. In this 
last regime, the spin down channel is also pinned at uj ~ 
£f = 0, and the spin up is displaced from it. Note that, 
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for this set of parameters, the Fano-Kondo dip splitting 
is resolved in the LDOS. The inset (II) displays more 
precisely the spin-polarized antiresonances analogous to 
those observed in Fig. 1(b) of Ref. [H|. In Fig. [S^b) we 
present the SR-LDOS as a function of the energy u> for 
different lateral STM-probe distances. As fcpi? increases, 
the Fano-Kondo dips disappear gradually, and for fc^i? ~ 
10.0, we recover the uncorrelated conduction band (see 
the inset (III)). Similar behavior is observed in Fig. (3) 
of Ref. El. 



VI. CONCLUSIONS 

We studied the spin-resolved local density of states 
(SR-LDOS) of a system composed of a Kondo adatom on 
a spin-polarized two dimensional electron gas in the pres- 
ence of an unpolarized STM probe. We derived the SR- 
LDOS expression in both the presence and absence of the 
STM probe. Our expression to the SR-LDOS in terms 
of phase shifts is general and independent of the method 
employed to calculate the local adatom GF. To deter- 
mine this GF we used the atomic approach in the limit 
of infinite U. The coupling parameter between adatom 
and host is assumed constant (local coupling). 

We were able to study the SR-LDOS in all the inter- 
ference regimes, varying the Fano factor q a . The main 
effect of the polarization was the tendency of one spin 
peak (90 = 100) or dip (go = 0.01, qo = 1) in the SR- 
LDOS to remain pinned around the host Fermi level as 
the polarization P increases. In contrast, the other spin 
peak or dip is shifted and loses amplitude as P increases. 
This contrasts to the usual behavior in the presence of 
a magnetic field, where the Kondo resonance is symmet- 
rically spin split and destroyed as the magnetic field in- 
creases, while here it is enhanced in one channel (down) 
and destroyed in the other (up). 

Our simulations are in close agreement with recent 
experimental results on adatom coupled to a ferromag- 
netic island.— The present system is a potential candi- 
date to promote the Kondo peak splitting without ap- 
plication of huge magnetic fields, necessary for adatom 
systems characterized by a large Tk- In particular, for 
the Fano factor q a = 0.01, we observed a continuous sec- 
ond order insulator-metal transition driven by the po- 
larization as presented in Fig. Ela). Finally, our model 
was able to describe qualitatively several experimental 
results^^L^ 
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Appendix A: The atomic GF for the Kondo adatom 

In this Appendix we present expressions employed in 
the spin dependent GF for the Kondo adatom considering 
the SIAM in the atomic version given by Eq. (|13|) of 
Sec. iLIBl To obtain the atomic GF, we use Zubarev's 
notation^ 



< f I Xg.^ \j > 
v-iEi-Ey) 



(Al) 



where /3 = l/fc^T with as the Boltzmann constant, 
T is the system temperature and O is the thermodynam- 
ical potential. The eigenvalues Ej and eigenvectors | j > 
correspond to the complete solution of the SIAM rlamil- 
tonian. The final result is the following 



uj - m 



(A2) 



where the poles and the residues are 

u\ a = E 3a — E la = E Sa — E 5cr = E-ic — E ia 

= ~ {Sdka — 5 a ) 
«2cr = E 5rT — E\ a = E$ a — Es a = E-! a — E 2a 
— - (Sdka + 8 a ) 
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(A4) 



respectively, which are defined in terms of 



(J 


— Ed — ef 




= £ko — £f 


Sla 


= sin 


C\a 


= cos 4> a 


&2a 


= sinA CT 


C2c 


= cosAcr 


(,1 


+ Eka — Sdk 



6a = ((e k „-e d ) 2 +4V 2 ) 1 / 2 
S'^de^-edf + SV 2 ) 1 / 2 

with ef being the Fermi energy and 

2V 



(A5) 



tan^ = 
tanA ff = 



Eka — £d + 5a 
2V2V 



Eka - £d + S' a 

Appendix B: Fano factor for the FM host 



(A6) 



In order to determine the Fano factor given by Eq. (|3"0)l 
in Sec. IIII Al due to the adatom-host coupling, we extend 
the procedure proposed for unpolarized bulk electrons^ 2 , 
to conduction states of a FM surface. To that end we per- 
form the calculation assuming the wide-band limit con- 
ditions lj <C D a and T <C D a in the advanced GF 



Go (w, R) 



Jk.R 



N F Ma ^P + E^W- E ka ~ if] ' 
k 



(Bl) 



with 77 — ^ + , which allows, in combination with Eq. 
(|3"2")l . to establish the following equalities 

5R{G CT (w,i?)} = U{g a (uj,R)} (B2) 

and 

9 [G a («, R)) = -3 {g a (u, R)} = TT Po A a (R) . (B3) 
Considering 

1 f 2v 

J (kR) = — exp[ikRcos9 k ]d6 k (B4) 
2?r Jo 

as the angular representation for the zeroth-order Bessel 
function in Eq. (|32p. and according to Eqs. (f3T)|) and 
(IB2D. the Fano factor becomes 



(B5) 



We can calculate this parameter rewriting the previous 
equation as 



QFMa = G CT (w, R) - iA a (R) , (B6) 



noting that the amplitude A a (i?) is already known from 
Eq. (jl2l). 

Thus, the quantity -^G a {w, R) must be found to pro- 
vide the relationship for the Fano parameter, which can 
be done using the decomposition 



1 ~ i jjy 1 PFMc 

npo 2 p 



written in terms of the integral 



Y,G,(u,R) (B7) 



1=1 



G t (u,R) = i / dskaHh 
1 



(0 



"F<T 



£fec 



R 



r 2 



Zko - iri 



(B8) 



that depends on the Hankel functions Hq\z) = Jo{z) + 
iY (z) and H '(z) = Jq{z) — iY (z). We conclude that 
the task here reduces to evaluate the integrals G\{lo,R) 
and G 2 (lj, R). 

The first integral is calculated by choosing a counter- 
clockwise contour over a semi-circle in the upper half 
of the complex plane that considers the simple pole 
e^ = +«r. Following the residue theorem, we have 



G 1 (w,R)=Hb 



(i) 



k Fa [ 1 + 



r 



lj — iT 



(B9) 



For the evaluation of the second integral we used a clock- 
wise contour over a semi-circle in the lower half plane 
with poles placed at Eka- — lj — if] and Eka = —^, which 
leads to 



G 2 (w,R)=2iH% 



(2) 



xH, 



(2) 



k F a 




-) 


R 


T 2 




r 2 + uj 2 


k F a 




r 1 

i — 

Da, 


)* 





lj + iT 
(BIO) 



As the complex conjugate property -Hq (z*) = 
(2) * 

Hq (z) is valid we are able to derive Eq. ([317)) from 



Eq. (HH considering Eqs. (|B7|) . (|B9]l and (|B10j) in the 
wide-band limit characterized by the conditions lj <C -Dcr 
and T < D a . 
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